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TECHNICAL NOTE 3508 


LAMINAR FREE CONVECTION ON A VERTICAL PLAIE WITH ERESCEIBED NONUmFORM 
NALL HEAT FLUX OR ERESCEIBED NONUNIFORM WALL TEMPERATURE 

By E. M. Speirrow 


SUMMARY 

An analysis is made for laminar free convection on a vertical plate 
with nonuniform thenaal conditions at the surface. Prescrihed variations 
are considered for the wall heat flux and for the wall temperature. 

For the situation where the wall -heat -flux variation is prescribed, 
graphs are presented from which the resulting wall -temperature variation 
may he obtained. Local heat-transfer coefficients may he readily deter- 
mined using the information given on the graphs. Results for the impor- 
tant special case of uniform wall heat flux are also given. 

For the situation where the wal 1-tenperature variation is prescribed, 
graphs axe presented from which the over-all rate of heat transfer from 
any length of the plate may he obtained. Another set of graphs is pre- 
sented for obtaining local heat-transfer coefficients . 

All the aforementioned results are given for fluids having Prandtl 
numbers in the range 0.01 to 1000. 

The flow is taEen to be of the boundary-layer type, and the problem 
is formulated by the Karman-Pohlhausen method. The solution of the resiilt- 
ing equations is achieved by series expansion. The first term of the 
series corresponds to the result for uniform thermal conditions on the 
wall. The succeeding terms give the influence of the nonuniform thermal 
conditions. The first five terms of the series have been calculated. 


INTRODUCTION 

T.«mi-naT- free convection on a vertical plate has been a subject of 
study since 1881. Most of the analytical work has been done for the sit- 
uation where the wall tenperature is uniform over the entire surface. 

An exact solution of the boundary-layer differential equations for free 
convection on a vertical flat plate with uniform wall tenperature is 
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given in, reference 1 for several Prandtl numbers in the range 0.01 to 
1000. Reference 1 also lists other \rork. on the unif orm.-'wal 1 -temperature 
problem. The K^rman-Pohlhausen method is applied to the uniform-wall- 
ten^erature case in references 2 (pp. 671-673) and 3 (pp. 158-163). The 
results reported in these references agree well with those from the ex- 
act solution of reference 1. 

Considerably less analytical wort has been done for the situation 
where the heat transfer is uniform over the surface. An exact solution 
of the boundary-layer differential equations for the uniform-heat -flux 
case is given in reference 4=, which also contains results calculated by 
the Kaiman-Pohlhaiisen method. The results of the exact solution and 
those from the TC»-rmRTi -Pn b1 bw isen method are in good agreement. 

Accounts of experimental investigations of free convection on ver- 
tical surfaces are given in references 1 to 6. 

In a large number of technical applications the ther mal conditions 
on the s\rrface are nonuniform. ' These nonuniformities in thermal condi- 
tions may be grouped into two categories: 

(1) The heat flux; may be prescribed to vary over the surface. It 
is then of interest to calculate the resulting variation of the surface 
tenqperatrxre . 

(2) The variation of the temperature on the surface may be pre- 
scribed. It is then of interest to calculate either the local rate of 
heat transfer at various locations on the surface^ or the over-all rate 
of heat transfer from the surfane, or both. 

This report presents a first attempt at solution of the free- 
convection problem on a flat plate for these two categories of non- 
viniform thermal conditions at the surface. 

This analysis was made at the HACA Lewis laboratory. 


gRTTCRAT. CONSIDERATIOHB 

The physical problem and the coordinate system are indicated in the 
following sketches, -vdiich show a vertical surface that may represent a 
flat plate or a vertical cylinder of large diameter: 
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Sketch (a). - Heat transfer from Sketch (h). - Heat transfer frcm 
wall to flTild and > t fluid to wall and t < t 

everyv;here . everywhere ^ 

Sketch (a) depicts a sitiiatlon in which the heat transfer i^ from 
the ^/all to the fluid at all points on the surface and the local wall 
ten 5 )erature t^ everyvdiere exceeds the aiabient temperature t^^. On the 
\tbAJ., either the heat flux may he designated to vary with x, or the 
wall ten^jerature may he prescribed to vary with x. (No variations are 
considered in the direction of the z-coordinate, normal to the pagej 
the problem is thios taken to he two dimensional.) The fluid in the 
neighborhood of the wall has a hi^er teu 5 )erature and a lower density^ 
than the flTiid far from the wall. Thi:iB, because of buoyancy, there will 
be established an up^'rard flow of fluid in the neighborhood of the wall. 
The region of space in which the upward flow primarily occurs is called 
the velocity boundary layer. A thermal boundary layer is defined as that 
region of space where the ten^ierature t deviates markedly from the am- 
bient ten^ierature t^^. In general, the velocity and thermal boundary 

layers have different thicknesses, the relative magnitudes depending upon 
the fluid properties . Both boundary layers are assumed to have zero 
thickness at the leading edge (x = O) . The velocity boundary layer is 
sho™ schematically in the sketch. 

Sketch (b) shows a situation in which the heat transfer is from the 
fluid to the tTall at all points of the surface and tlie local wall ten^ier- 
ature t^ is eveiywhere less than t^. Again, on the wall, either the 
heat flijx or the ten^jeratiire may be prescribed to vary with x. Here, 
the flow of fluid in the boundary layer is downward as shown. 


■^Dhis refers to fluids showing the iisual trend of density decreasing 
trLth increasing temperature. 
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If the coordinate systems are taken as shown in the sketches^ the 
method of analysis and the results for the heat -transfer parameters are 
the same for these two situations; and there will he no need to treat 
them separately. So, the analysis will he carried out for the case of 
heat transfer from the wall to the fluid (t^ > t^^), hut it is to he remem- 
bered that the results apply to hoth situations depicted in the sketches. 

A definite class of variations of the wa l l heat flux will he pre- 
scribed'. Suppose that the region of interest on the vertical plate lies 
between x = 0 and another location x = x^^ (xj^ must lie in the region 

of laminar flow over the plate). The VTall heat fluxes to he considered 
here have a finite nonzero value at x = 0 and either increase steadily 
from X = 0 to x = Xj^ or else decrease steadily from x = 0 to 

X = Xl. 

The form of the wall -temperature variations considered here is sim- 
ilar to that outlined for the wall -heat -flux variations. The wall tem- 
perature relative to ambient has a finite nonzero value at x = 0 and 
either rises steadily frcm x = 0 to x = Xj^ or else decreases stead- 
ily from X = 0 to x = xj^. 

Although the analysis made here is for free convection in a gravity 
fields it may easily he generalized to include other force fields. It 
js only necessary to replane the gravitational force per unit mass g 
by the body force per unit mass of the other force field imder 
cons iderat ion . 


BASIC EQUATIONS 


The equations expressing conseivation of mass, mcanentum, and energy 
for steady laminar flow in a boundary layer on a vertical flat plate axe 
eis follows: 


Su dv _ 


^ du „/. 

u ^ + V — = gp(t 

ox oy 


tj + V 


u T h V — = a 


dx 




dy^ 



( 1 ) 

o 

( 2 ) 


( 3 ) 


(ah symbols are defined in appendix A.) In accordance with the usual 
practice in free convection, the density is considered a variable in 
formulating the buoyancy term gp(t - tg^) . Aside from this, the fluid 

properties are taken constant. Viscoias dissipation and work against 
the gravity field are neglected. 
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Following reference 3, it is assimed that a common boundary-layer 
thdckness 6 can be used for both velocity and thermal bovindary layers. 
This assumption has its justification in the fact that the resiilts of 
calculations performed with it are In good agreement ^d-th those from 
exact solutions of the boundary-layer differential equations for the 
cases of unifoim wall ten^erature and uniform heat flux. 


Then^ eqmtions (l) to (3) are integrated across the botuadary 
layer to give 



_d^ 

dx 



u(t 



(3a) 


The Integrated form of equation (l) has been absorbed into equations (2a) 
and (3a) . 

These equations have a definite physical meaning. They are, in fact, 
expressions of the conservation laws for the element of boundary layer 
shown in the foUovring sketch: 



Sketch (c) 

Solutions for the velocity and tengjerature in the boundary layer 
will be obtained Tdiich satisfy these conservation equations and the 
boundary conditions. These solutions will in turn be used to calculate 
the imp ortant heat -transfer parameters. 
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ANALYSIS 


Prescribed Nonunifo3rm Wall Heat Flux 


The analysis proceeds with the use of the Karman-Pohl ha\isen method, 
according to -vjhich the velocity and teiqoerature distributions in the 
boundary layer are written an polynomials in y whose coefficients are 
functions of x. The coefficients are found from the boundary conditions 
of the problem and by using the Integrated momentum and energy equations 
(eqs. (2a) and (3a)). 


The following polynomials are chosen: 


u 




U) 

(5) 


&t 


Equation (4) satisfies the conditions: when y = 0; t = t^^ 

and -^ = 0 (smooth-fit condition) when y = B. The wall heat flux q 
is to be regarded as a specified function of x. The conditions satis- 
fied by eqxxation (s) are: u = 0 when y = 0} u = 0 and = 0 

(smooth-fit condition) ^dien y = 6. The functions ai(x) and B(x) still 
remain to be determined. 

The polynomials representing the velocity and ten 5 )erature distribu- 
tions are introduced into equations (2a) and (3a), and after the inte- 
gration is carried out, there results a pair of first order, ordinary 
differential equations for o and 6. In (dimensionless foim, these 
eqmtions are 


105 dK 


(2"a) = %- 


2 

A 


( 6 ) 


30 dX 


^2 _a. 

oJ 


%J ^ % 


(7) 


X, 2, and A are the dimensionless counterparts of x, O), and 8 and 
are defined in the Symbol list (appendix A) . The symbol q^, which 
represents the heat flux at the leading edge (x = O) , is used as a 
reference heat flux in the rest of the analysis. 
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next, the variation of q/q^ ie to he prescribed. Suppose that 
the region of Interest on the vertical plate lies between x = 0 and 
another location x = x^. The surface heat fluxes to be considered here 
have finite nonzero values at x = 0 and either Increase steadily from 
X = 0 to X = Xj^ or else decreeise steadily from x = 0 to x = Xj^. 
Explicitly, the class of surface-heat -flux variations to be considered 
is -v/ritten as 


t = " " ‘(tf 

The exponent r, Tdiich gives the shape of the variation, is retired to 
be a positive number. Integral or nonintegral. From equation ^8)^ i't 
may be seen that e represents the maximum (percentage) deviation of 
q from q^ in the region x = 0 to x =? Xj^. 

For the q/q^ given by equation (s), the pair of differential 
equations (S) and (?) can be solved by expanding 2 and A in Maclaurin 
series in terms of e as follows: 

2(X,6,Pr,r) = ^ e^^(X,Pr,r) = 2Q(X,Pr,r) 

n=0 

00 

A(X,6,Pr,r) = ^ e“^(X,Pr,r) = ^(X,Fr,r) 
n=0 

The expressions for 2, A and q/q^ are Introduced into equations 

(e) and (7). Terms are gro\;q)ed according to the power of s that mul- 
tiplies them, that is, 

]■" ■■■ 

In order that equations (ll) and (l2) be satisfied for any value of 6, 
each of the brackets must' be identically zero. Equating to zero the 
brackets multiplying s® yields a pair of simultaneous equations for 
2 q and Aq, In a similar fashion, the brackets multiplying give 

a pair of equations for 2^ and A^. It may be noted that the equations 
for 2^ and A^ will include all the 2 q through and the A^ 

through A^_ 2 _ • 


+ e2jL(X,Pr,r) + ••• (9) 
+ e-^^(X,Pr,r) + •••(lO) 
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The solution for 2 q and which coincides with the results 

for uniform wall heat flixx (e = O), is treated more fully in appendix 
B. The succeeding 2^^ and evidently give the effects of the non- 

■uniformity of the heat flTox. 

The first five functions in the Maclaurin series (2g through 2^ 
and Aq throu^ A^) have been con^juted. The results are listed in 
concise form as follows: 


= (±1)“ X^(Pr,r)^5 


(13) 

= (±l)° rJ.Pr,x)fl 

L\“‘x^A 

(14) 


The factors and were calculated from linear algebraic equations. 

The results thus found are applied in a later section in which the 
heat -transfer parameters are obtained. 


Prescribed NonuniCorm Wall Ten^jerature 

The analysis for the prescribed-nonuniform-wall- temperature case 
proceeds in a similar fashion to that for the prescribed heat-flux case. 
The te^perat^^re and velocity profiles in the boundary layer are approxi- 
mated by the following polynomials; 

t - ta= ' i) ® - i) 

u . ® I (l - 2) (16) 

Equation (l5) satisfies the conditions: t = t when J = 0} t = t„ 

^ t- 

and -^ = 0 when y = 6. The wall tenperature t^ is to be regarded 

as a specified function of x. Since the ambient tenperature tg^ is 
talen to be constant, then 0 is regarded as a prescribed function of 
X. The conditions satisfied by equation (l6) are: u = 0 when y = 0; 

u = 0 and ^ = 0 when y = 6. 

The variation of the tjall tenperature t^ is prescribed in a re- 
gion of interest between x = 0 and some other locauion x = x^. The 

wall tenperature relative to anbient considered here has a finite non- 
zero value at x = 0 and will either increase steadily from x = 0 to 
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X = Xj^ or else decrease steadily from x = 0 to x = Xj^. The explicit 
form of the tenperature variation considered here is 



Again,, the exponent r is a positive number, integral or nonintegral. 

As may be seen from equation (l7a), c is the maximum (percenteige) de- 
viation of 9 from 0 q in the region 0 < x < x^. 

Equations (l5) and (l6) are introduced into the Integrated momentum 
and energy eqmtions (eqs. (2a) and (3a)). The resulting simultaneous 
equations for OJ and 8 are solved subject to the prescribed wall- 
teraperature variation (eq. (l7)) using Maclaurin series. The first five 
functions in the Maclavirin series have been calculated. The results 
have a form similar to those given in equations (l3) and (l4) . The 
solution for the uniform-wall -temperature case coincides with the lead- 
ing term of the Maclaurin series expansion. 


HEAT-TRARSFER KESUL3B 
Prescribed Nonmlform Wall Heat Flux 


For the case where the wall heat flux is prescribed, it is of in- 
terest to determine the resulting wall-temperature variation and local 
heat -transfer coefficients . 


Res\ilting t'fa.n -temperature variation . - The wall tenperature is ob- 
tained by setting y equal to 0 in equation (i). Thus, 


t - t 

V a - 2k 


(18) 


or 


w 




(18a) 


For a imiform wall heat flux q = q^, it has already been noted that 
A = dQj so. 


' w 


^a)q 


2k 


(19) 
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■v^iiere (t - t ) is the ■wall-teraperature variation corresponding to a 
w a- 

uniform heat flux q = <1^. In equation (l8a), is evaluated from 
equations (8), (lO), and (l3), and q^5o/2k is evalviated from equation 
(l9), giving the result 



( 20 ) 


2 

OJhis equation gives the ratio of the wall teirperature at some loca- 
tion X (x < having a specified ratio q/q^ to the wall temperature 

at the same location on a plate having a uniform heat flux q = q^. 

Since the X's depend upon the Prandtl number and r, the temperature 
ratio given hy equation ^20 ) depends on q/q^^ Pi* and r. 


The ratio 


(t. 


w 




^^w " ^a^qo 


is plotted in figures l(a) to (e) . Each of 


the plots applies for a specific value of r. The values of r = 0, l/2, 
1^ 2f and 3 have heen used for the five plots. Results for other values 
of r may he obtained by replottlng^ the information given in figures 
l(a) to (e) using r as the abscissa variable. On each of the plots, 
the tenperature ratio is plotted against q/q^ for the range 


0.5 < q/<l^< 1.5 with Prandtl numbers between 0.01 and 1000 appearing 
as parameters on the c;xrves. In cases where curves for different Prandtl 
numbers fall so close together as to make it inpossible to plot them 
separately, one curve was -used for the several Prandtl numbers. 


Once the ratio /. - — qr-f — has been determined from one of the 

figures (or by replotting the data given therein), the wall temperature 
(t^ - tg^) can be found if an expression for (t^ - ^a^q^ given. The 


expression for (t^ - ^a^q^ 

case of uniform heat flux 
Karmlm-Pohlliausen method. 


, the wall temperature correspond i ng to the 

q = q^, is derived in appendix B by the 
and appears on figure 1. 


^'Relative to ambient" is to be tinderstood every time the wall tem- 
peratirre is mentioned. 

^r = 0 is a limiting case, which is included here to facilitate 
the replottlng when results for r near zero are required. 


3430 
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Local heat -transfer coefficient . - The definition of the local 
heat -transfer coefficient is 


OJ 

(F- 

CM 

o 



For the case of a uniform heat flux q = q^, equation (2l) is 

h '^O 

K - 

The ratio of equation (2l) to equation (21a) is 


( 21 ) 


(21a) 


( 22 ) 


This equation gives the ratio of the local coefficient at some position 
X having a specified ratio q/q^ to the local coefficient at the same 
position on a plate having a uniform wall heat flux q = q^. When o/q^j 

Pr and r are known, the ten^peratirre ratio appearing in equation (22) 

can be calculated. The expres- 


XX CLLJ.*.!. X CLX ^ rv i it. iw I LJ^Ci-L U LXX < 

is known (from fig. l)j hence, c 

Sion „ is given in appendix B. 


'Prescribed Nonuniform Wall Ten^jerature 


For the case where the wall teng)erat\are is specified, it is of in- 
terest to calculate the over-all heat-transfer rate and local heat- 
transfer coefficients . 


Over-all heat-transfer rate . - The over-all rate of heat transfer 
in a section of plate of width b frcm x = 0 to another location 
x(x <: Xj^) is 



(23) 


The derivative (^/^)y^ is evalmted from equation (l5) . 
into equation (23) gives 


Q = 2kb 


/ 


" ^a^ 


dx = 2kb 




dx 


Substitution 


6 


(24) 
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For the case of imiform wall temperatxire 0 = 6^, eqixatlon (24) hecomes 
The ratio of equation (24) to equation (25) is 



Equation (26) is the ratio of the over-all heat transfer from 
X = 0 to X (x < Xj^) on a plate with variable surface temperature to 
the heat trajasferred in the same region on a plate with a uniform wall 
ten^ierature 0 = .0^. Since 8 corresponding to the prescribed wall- 
ten^ierature variation (eq. (l7)) has been calculated in the Analysis 
section^ it is possible to ev8Ll\iate the integrals appearing in equation 
(26). 

The results for Q/Q^ , 'vjhich are shown in figures 2(a) to (e); 

depend upon the value of location x^ upon Prandtl num- 

ber,, and upon r. Each plot is for a specific value of r. The values 
of 0^ 1 / 2 , 1, 2, and 3 have been used for the five plots. Results for 
other values of r may be obtained by replotting the information given 
in figures 2(a) to (e) with r as the abscissa variable. On each plot, 
the ratio Q/Qq is plotted against for the range 0.7< 0/0 q< 1.3, 

with Prandtl numbers bet^reen 0.01 and 1000 appearing as a parameter on 
the curves . 

It is to be emphasized that the value of s-'t x is to be 

used for determining values from the figure when the over-all heat trans- 
fer from X = 0 to x is required. The quantity Qn , which is the 

o 

over-all heat -transfer rate frcm x = 0 to x on a plate having uni- 
form surface temperature 0 = 0^, is calculated in reference 3 by the 

Kaim^-Pohlhausen method. The expression for is given in 

figure 2. ° 

Local heat-transfer coefficient . - Prom the defining equation (eq. 
( 21 )), the local coefficient is found to be 
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The derivative is obtained from equation (l5) . For mlform wall teng)er- 
ature 0 = therefore 

where tu a is the local coefficient for the case of a iinlform wall 
temperature 0=0^. Combining equations (27) and (28) gives 

This latio is plotted in figures 3(a) to (e) in the manner already out- 
lined for the preceding results of the analysis . The expression for 
hjj. Q (from ref. 3) is given on the figure. 



DISCUSSION 


Ohe following generalizations can be made from inspection of the 
results ; 


(l) At a location x with a fixed the deviation of 

^'w “ ■ta 

-rr r— T — from 1 increases with decreasing r for a fixed Prandtl 

“ "^a^qo 

number. Also, the -deviation from 1 is larger for larger Prandtl numbers 
at a fixed r and q/^Q* 


with 
1 is 


(2) In a length of plate between x = 0 and x, having a fixed 
at X, it is seen that the deviation of Q/Qq from 1 increases 

decreasing r at a fixed Prandtl niimber. Also, the deviation from 
smaller for the larger Prandtl nimbers at a fixed r and 0 /®q- 


(3) The deviation of 

at a fixed Prandtl nimiber 
from 1 is smaller for the 




frcm 1 increases with increasing r 


and a fixed q/lQ at x. Also, the deviation 
larger Prandtl numbers at a fixed r and q/qQ* 


(4) The deviation of h^h^^ q from 1 increases with increasing r 

^ o 

at a fixed Prandtl number and a fixed at x. Also, the deviation 

from 1 is smaller for the larger Prandtl numbers at a fixed r and 

e/So- 
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It has already been noted that the tenrperatxire ratio 






■vrais plotted in figure 1 for 0.5 ^ q/q|^-^1.5. This range was decided 

upon by studying the extent of the errors due to tinncation of the 
infinite series in equation (20) after n = 4. The range 
0.7 < 0/0Q< 1.3 used in figures 2 and 3 was decided upon by studying 
the truncation errors of series associated vrith equations (26) and (29). 


COISCLUDIMG EEMAEKB 

When the nonuniform, wall heat fltax can be written in the form of 
equation (8), the resiilting wall teE5>eratures may be found directly from 
figure 1, and the local heat-transfer coefficient can be calculated from 
equation (22). All that is needed to use the graphs is the value of 
q/q^ at the point of interest^ the exponent r which gives the shape 
of the q/q^ variation, and the Prandtl number. It is to be noted that 
E is not needed. 

When the wall temperature is specified by a relation of the form of 
eqtiatlon (l7), the over-all heat-transfer rate and local coefficients are 
fomd from figures 2 and 3, respectively. The over-all heat -transfer 
rate for a section of plate from x = 0 to x (x Xj^) may be found from 

the graphs when the following are knoT'/n: ®/®o exponent r 

which gives the shape of the variation, and the Prandtl number. 

The same quantities are needed to find the local heat-transfer coefficient 
at X. 

It is recognized that whenever a new application of the Karman- 
Pohlhausen method is made, it is desirable to confirm the resiilts by 
checking with those of experiment or of a less approximate analysis. 

The author is not acquainted with any experimental data that may be 
used to check the results derived here for nontinlform thermal conditions 
at the siirface. Uor is there now available any other analysis with tiiich 
the present results may be ccmpared. 

For the special cases of uniform wall temperature and uniform heat 
fltix, there are exact solutions of the laminar-boundary-layer equations 
as well as experimental dateu For these cases, the heat-transfer re- 
sults derived from the Kerman -Pohlhausen method agree well with those 
frcm the exact solutions and those of experiment. 


Lewis Fli^t Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 4, 1955 
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APPENDIX A 


SXMBOIS 


The follOTring symbols axe used in this report: 


vri-dth of plate, ft 

specific heat at constant pressure, Btu/(lh)(‘^) 
acceleration due to gravity, ft/sec^ 

modified Grashof nuniber based on x, 5—, dimensionless 

3 k'’ 

X 

Grashof number beised on x, — , dimensionless 

local heat-transfer coefficient, Btu/(sec) (sq. ft)(‘^) 

thermal conductivity, Btu/(sec) (ft) (°P) 

index for naming terms of a series, dimensionless 
c Lt 

V 

Prandtl nxmiher, , dimensionless 

^ k a' 

over-all heat-transfer rate on a plate of vidth h Be- 
tween X = 0 and x (x < Xj^), Btu/sec 

loceil heat-transfer rate per unit area, Btu/(sec)(sq ft) 

heat-transfer rate per unit area at x = 0, Btu/(sec)(sq ft) 

exponent defined hy eqs. (8) and (l7), dimensionless 

static tenperature, 

velocity component in x-direction, ft/sec 

velocity conponent in y-directlon, ft/sec 

coordinate meeisuring distance along plate from leading 
edge, ft 

coordinate defining region of Interest x = 0 to x = Xj^, 

/ 

f gPqoV^^ 

dimensionless coordinate, I s x 
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y 

a 




■^0 


0 

0 


'0^ 


V- 

V 

o 

OJ 


■ T. 


^ • • A 


V 


2o^2jl^ * 


•2, 


coordinate measuring normal distance from plate, ft 


thermal diffusivity, , sq ft/sec 


PCt 


coefficient of thermal expansion, (°F) 


-1 


fantors in solutions for — 2^ defined by eq. 

(l4), dimensionless 

dimensionless houndary-layer thickness for prescrihed- 

1/4 


heat-fl\3x case 


■m 


coefficients in Maclaurin series expansion of A ±q 
terms of £, dimensionless 

boundary -layer thickness, ft 

houndary-layer thickness for uniform heat flux or uniform 
wall ten 5 )erature, ft 

nimiher giving percentage deviation of q at x = Xj^ from 
q^, or number giving percentage deviation of 0 at 
X = x^ from 0^, dimensionless 

wall- to ambient-tenperature difference, t^ - t^^, 

wall -to amblent-teiig)erature difference at x = 0, °F 


•A^ defined by eq. 


factors in solutions for , 

(l3), dimensionless 

absolute viscosity, lb/(sec)(ft) 

kinematic viscosity, sq ft/sec 

density, Ib/cu ft 

velocity fimction defined by eqs. (s) and (l6), ft/sec 

dimensionless velocity function, 1 g — j CD 

coefficients in Maclaurin series expansion of 2 in 
terms of e, dimensionless 
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Siibscripts; 

a ambient 


on a plate having uniform heat flux 
w wall 

6 ^ on a plate having uniform wal 1 - to ambient-teraperature 

difference 6 ^ 
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APPENDIX B 


EESUDIS FOR UNIFORM VZALL BEAT FLUX 
When for all values of x, equations (6) and (?) become 


JL A (aV ) . ^ ?o 

105 dX “6 


1 d- /q a^\ _ 2 


(Bl) 


(B2) 


The 


solutions for Sq and Aq are 


2q = ( 6000 )^/^ Pr“^/^^| + Pr^ ^ X^/' 


Aq = (360)^/^(| + Pr) Pr“^/^ X^/^ 
Equation (B4) can be rewritten as 

,1/5 

^ ^ - (360)V5| 


X 


(I * 


J 


where Gr* is a modified Grashof number based on x and defined by 


(B3) 

(B4) 

(B5) 


Gr* = 

X 




kv 


(B6) 


The surface -ten 5 )erature distribution found by Introducing equation (B5) 
into eqmtion (l8) is 


(t, - t )„ = 1.622 

^ w a'q^ 


f O.8 + Pr\ ^^^ 
^ yPr^ Gr* j 


(B7) 


The local heat-transfer coefficient is obtained by the foUo^rtng 
rearrangement of eqiiation (B7): ^/_ 

qp 0 . 62 k f P^^g^x' \ 

^^^o (-^w - -^a^qo ^ 


(B8) 


I 
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Figure 1. - Ratio of xall tempera txzre at acme location x on plate with variable heat flux to vail 
teiDperature at same location on plate with imlform beat flux q « . Abscissa is value of 

<i/^ at point of interest} r Is an exponent In equation (e)(relatioa by which heat-flux varla- 
tl^ is prescribed); Prandtl number is parameter on curves. 
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FlSuro 1, - Continued. Ratio of wall tampcrature at sono location x on plate with varlatle heat 
flujc to wall tooporaturo at ea*» location on plate with unifora heat flui q " ' AhsoloBa lo 

value of q/On at point of Interant; r 1b an exponent In equation (5) (relation oy Mhloh heat- 
flux varlatlofl 1b pireeoribed) ; Prandtl maiiber IB parameter on ourveB. 
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(q) P - 3. 

t«np«rature at sooe location x do plat* with verlabla heat 
fl« to wall te«ip*ratiiro at sana location on plate with uniform heat flux q - o_ . Abaolasa 1* 

exponent In equation (8) (relatlo^j whloh heat- 
flux variation la pr«Bcribe<l) j Prandtl maiibar la paranet^r on ourv«Be 
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(b) r - 1/2. 

Figure 2. - Continued. Ratio of over-all heat transfer from 
X = 0 to X on plate with variable wall tetpperature to over- 
all heat transfer In same region on plate v;lth uniform viall 
temperature 9 ■= 6^ . Abscissa Is value of 9 / 9 q at x; r Is 
an exponent Is equation (17) (relation by v;hlch wall-tempera- 
ture variation Is prescribed); Prandtl number Is parameter on 
curves . 


CZ-4 bact 
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Figure 2. - Continued. Ratio of over-all heat trans- 
fer from x = 0 to x on plate with variable wall 
temperature to over-all heat transfer In same re- 
gion on plate with uniform wall temperatTU?e 6 = 6 ^ 
Abscissa Is value of ^ exponent 

In equation ( 17) (relation by which wall -temperature 
variation Is prescribed); Prandtl number Is param- 
eter on curves . 
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Figure 2. - Concluded. Ratio of over-all heat trans- 
fer from X = 0 to x on plate with variable wall' 
temperature to over-all heat transfer in same re- 
gion on plate with unlfcwro wall temperature (9 t= 0 ) 
Abscissa is va l\^_ of S/Rq at x; r Is an exponent 
in equation (17 relation by which wall -temperature 
variation is prescribed); Prandtl number is param- 
eter on curves . 



30 


M.CA TN 3508 



(a) r « 0. 


Figure 3. - Ratio of local heat-transfer coefficient at x for 
plate with variable iiall temperature to local coefficient at 
same location on plate vjith uniform mil temperature 6 «» 6 
Jlbscissa is value of e/e^ at x; r is an exponent in equa- 
tion (17) (relation by which wall-temperature variation Is pre- 
scribed); Prandtl number is parameter on curves. 
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(b) r - 1/2. 

Figure 3. - Continued. Ratio of local heat-tr^sfer coefficient 
at X for plate with variable wall temperature to local coeffl 
dent at same location on plate with uniform vjall temperature 
„ « ni j 0/0 Tr. T. Is an exponent 


0 ^ 0 .. Abscissa Is value of >// Vq dU A.J X 
in equation (17) (relation by which wall-temperature variation 
Is prescribed); Prandtl number is parameter on curves. 
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Figure 3. - Continued. Ratio of local heat-transfer coefficient 
at X for plate vfith variable wall temperature to local coeffi- 
cient at same location on plate with uniform wall temperature 
0 = Oq • Abscissa is value of 0/0- at x; r is an exponent 
in equation ( 17) (relation by which wall -temperature variation 
is prescribed) ; Prandti number is parameter on curves . 
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Figure 3. - Concluded. Ratio of local heat-transfer coefficient 
at X for plate with variable wall temperature to local coeffi- 
cient at same locatlcsi on plate with uniform wall temperature 
e - . Abscissa Is value of 0/0o at; x; r Is an exponent 

In equation (l7)(relatlon by which wtai-temperature variation 
Is prescribed); Prandtl number Is parameter on curves. 
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